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Abstract 

In scenarios with extra dimensions and TeV-scale quantum gravity, black holes are ex- 
pected to be produced in the collision of light particles at center-of-mass energies above the 
fundamental Planck scale with small impact parameters. Black hole production and evapo- 
ration may thus be studied in detail at the Large Hadron Collider (LHC). But even before the 
LHC starts operating, neutrino telescopes such as AMANDA/IceCube, ANTARES, Baikal, 
and RICE have an opportunity to search for black hole signatures. Black hole production 
in the scattering of ultrahigh energy cosmic neutrinos on nucleons in the ice or water may 
initiate cascades and through-going muons with distinct characteristics above the Standard 
Model rate. In this Letter, we investigate the sensitivity of neutrino telescopes to black 
hole production and compare it to the one expected at the Pierre Auger Observatory, an air 
shower array currently under construction, and at the LHC. We find that, already with the 
currently available data, AMANDA and RICE should be able to place sensible constraints 
in black hole production parameter space, which are competitive with the present ones from 
the air shower facilities Fly's Eye and AGASA. In the optimistic case that a ultrahigh energy 
cosmic neutrino flux significantly higher than the one expected from cosmic ray interactions 
with the cosmic microwave background radiation is realized in nature, one even has discovery 
potential for black holes at neutrino telescopes beyond the reach of LHC. 



1. The possibility that black holes will be produced in the collision of two light particles at 
center-of-mass (cm) energies above the Planck scale with small impact parameters PP seems now 
within reach in the context of theories with 6 = D — A > 1 large compact 2, or warped [3] extra 
dimensions and a low fundamental Planck scale Md > 1 TeV characterizing quantum gravity. In 
these theories one might expect the copious production of black holes in high energy collisions at 
cm energies above Md e4|. Correspondingly, the Large Hadron Collider (LHC) [Zj, expected to 
have a first physics run in 2006, may turn into a factory of black holes at which their production 
and evaporation may be studied in detail [SI El GDI HI] • 

Black hole production and subsequent decay in the scattering of ultrahigh energy cosmic neutrinos 
on nucleons in the atmosphere may initiate quasi-horizontal air showers far above the Standard 
Model rate ^12, 131 with distinct characteristics J3] . The reach of the Pierre Auger Observa- 
tory for extensive air showers, expected to be completed by the end of 2003, in the black hole 
production parameter range has been thoroughly investigated in Refs. [121 CHI EI • It was found 
that, depending on the size of the flux of ultrahigh energy cosmic neutrinos, the first signs of black 
hole production may be observed at Auger already before the start of the LHC. Furthermore, 
it was shown that already now sensible constraints on black hole production can be obtained 
from the non-observation of horizontal showers fHl ^] by the Fly's Eye collaboration [SHj and 
the Akeno Giant Air Shower Array (AGASA) collaboration pT], respectively. These constraints 
turn out to be competitive with other currently available constraints on TeV-scale gravity which 
are mainly based on interactions associated with Kaluza-Klein gravitons, according to which a 
fundamental Planck scale as low as Md = 0(1) TeV is still allowed for 5 > 4 flat or 5 > 1 warped 
extra dimensions 22^ 

In this Letter, we want to consider the sensitivity of neutrino telescopes such as AMANDA/Ice- 
Cube [23, ANTARES [21], Baikal [23, NESTOR [2E], and RICE [27] to black hole production 
in the scattering of ultrahigh energy cosmic neutrinos on nucleons in ice or water, considerably 
extending a first exploratory study [2H]- After briefly reviewing the phenomenological model 
of black hole production adopted, we calculate the rate expectation for events contained in the 
fiducial volume of a neutrino telescope. We then extend our calculations to rates from through- 
going muons. Finally, we discuss the achievable sensitivity for detection of black hole production 
and compare it with that from air-shower facilities and the LHC. 

2. Our phenomenological model of black hole production and decay will be based on the following 
working hypothesis which summarizes the current understanding of the production and decay of 
black holes in TeV-scale gravity scenarios [H El IHl El El UHl lUj ■ 

It is assumed that at trans-Planckian parton-parton cm energies squared, s ^ M|,, and at small 
parton-parton impact parameters, b <^ rs (Mbh = v^), i-e. at impact parameters much smaller 
than the Schwarzschild radius rs of a (4 + 5)-dimensional black hole with mass Mbh = |29j^ 



^For an exhaustive list of references in this context, see also Ref. 

^We define Md as in Ref. [Slj- Equation |^ is valid as long as rs <C Rc, with Rc being the compactification 
or curvature radii in the flat or warped scenario, respectively. 
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a black hole forms with a cross section 



a(^J _ bh) = a^^is) ^ vr r| (Mbh = 71) (71 - M^T) • (2) 

Here, M,^™ ^ M/) parametrizes the cm energy above which the semiclassical reasoning mentioned 
above is assumed to be valid. The contribution of black hole production to the neutrino-nucleon 
cross section is then obtained as the convolution 

^']^(^) = E rdx/,(a;,/.)CM. (3) 





Here, s = 2m]\fEy, with m^r = {rrip + m„)/2 denoting an isoscalar nucleon mass and the 
neutrino energy in the laboratory frame, denotes the neutrino-nucleon cm energy squared. The 
sum extends over all partons in the nucleon, with parton distribution functions fi{x,fi) and 
factorization scale fi. As in Ref. |16j, we have used in the calculation of the cross section Q 
various sets of parton distributions as they are implemented in the parton distribution library 
PDFLIB j^. In the following we display only the results obtained with the CTEQ3D 
parton distributions with /i = mm{\/§,10 TeV), keeping in mind that uncertainties associated 
with different parton distribution sets are in the (9(20) % range, whereas a choice of the other 
natural factorization scale /i = r^^ 01131 decreases the predicted rates by a factor of 0{2). 

It is fair to say that there remain nevertheless substantial uncertainties in the size of the cross 
section. As noted in Ref. [22] , the semiclassical production of black holes resembles largely the 
problem of baryon and lepton number violating processes ( "sphaleron [21] production" ) in multi- 
TeV ^ mw/ctw) particle collisions in the standard electroweak theory jHHl IHHl IHZI and 
more generally the problem of multi-particle production in weakly coupled theories [2H1, which is 
not yet completely understood [2ni- There might be an additional exponential suppression factor 
rendering semiclassical sphaleron, multi-particle or black hole production unobservable in the TeV 
range [221 i^^^y however, Ref. [101 )• Nevertheless, in this Letter we shall confine ourselves to the 
simple geometric cross section estimate Q and postpone an inclusion of a possible exponential 
suppression factor to a longer write-up [Hlj . 

After production, the black holes decay primarily via Hawking radiation , corresponding to a 
temperature 

into a large number [^, 

(n) ^ ^ — ^ , (5) 

of hard quanta (cf. Fig. GJ, with energies approaching several hundreds of GeV in the cm system. 
We shall assume that the average multiplicity ^ refers to Standard Model particles [SI - an 
assumption which is still under debate [121 ^ in a "flavor" -democratic fashion. 

3. Let us consider now the rate of black hole events taking place within the fiducial volume of 
a neutrino telescope. The number of black hole events initiated by neutrino-nucleon scattering 



3 



120 
100 
80 

S 60 
40 
20 


10^ 109 iQio iQii 10»2 
K [GeV] 

Figure 1: The average number (n) of Standard Model particles (0), into which a black hole 
produced in neutrino-nucleon scattering decays, as a function of the incident neutrino energy in 
the laboratory frame, E^. 




occurring per unit time t and solid angle with energy larger than a threshold energy Eth inside 
a subsurface detector volume V is (see e.g. Ref. |36j ) 

d'N^L pV r 



poo p , 

/ dE,F,{E,)al''^iE,) exp -al%iE,)X{e)/mN , (6) 



where p is the density of the material in which the neutrino interaction occurs, Fi, = +Fp.) 
is the sum of the differential diffuse neutrino fluxes, denotes the total interaction cross 
section as a sum of the usual Standard Model charged current cross section c'^j^ and the black 
hole production cross section (0), 

tot cc I bh ffj\ 

and, finally, X{9) is the column density of material between the detector and the upper atmo- 
sphere, as a function of the zenith angle 9. It can be approximated by p 



x{e) = p 



(i?e - cos2 9 + 2DR^ - - (i?® - D) cos 9 



where i?® denotes the radius of the earth and D is the vertical depth of the detector. 

As in Ref. |16.], we shall exploit both conservative lower and upper limits 1^1^ on the presently 
unknown differential flux F^y of ultrahigh energy cosmic neutrinos (for recent reviews, see Ref. |45p 
entering into event rates at neutrino telescopes such as Eq. (0). More or less guaranteed, and 
therefore comprising a reasonable lower bound, are the so-called cosmogenic neutrinos which are 
produced when ultrahigh energy cosmic protons inelastically scatter off the cosmic microwave 
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Figure 2: Predictions of the cosmogenic neutrino flux, F^, = J2i [Fui + Fp.]. Long-dashed (long- 
dashed-dotted) hne: Flux from Ref. jlH] for cosmological evolution parameters m = 2, z^aax = 2 
(m = 4, z^s.x = 4). Solid (dotted) line: Flux from Ref. jlH], assuming a maximum energy of 
^max = 3 ■ 1020(21) gV for the ultrahigh energy cosmic rays. Shaded: Theoretical upper limit 
of the neutrino flux from "hidden" astrophysical sources that are non-transparent to ultrahigh 
energy nucleons |44j . 

background radiation in processes such as J97 — A ^ nvr"'", where the produced pion subse- 
quently decays ^4Ij. Recent estimates of these fluxes can be found in Refs. jlHl HHl EOl EI] , some 
of which are shown in Fig. El The upper limit from "hidden" hadronic astrophysical source^ i. e. 
from those sources from which only photons and neutrinos can escape, is much larger jH] and 
also shown in Fig. |21 

The projected reach in the black hole production parameter space for 5 = 6 extra dimensions 
for contained events in an under-ice detector corresponding to the IceCube proposal |^ (2 km 
depth, flducial volume 1 km^) is shown in Fig. |H1 for two of the discussed ultrahigh energy cosmic 
neutrino fluxes. Table Q summarizes the corresponding expected Standard Model background 
rates for contained events due to charged current interactions. Taking into account the small 
effective volume, V ~ 0.001 -i- 0.01 km^, of the currently operating AMANDA and Baikal neutrino 
telescopes, it seems that, on the basis of contained events, these facilities cannot compete with 
the already established limits from air shower observatories such as Fly's Eye and AGASA, and 
the projected ones from Auger [121 CEl El • Even IceCube does not seem to be really competitive, 
as far as the contained events are concerned: According to the time schedule of IceCube, the 
construction starts in 2004 and the flnal effective volume, ^ 1 km'^, is reached only after the 

■^For an early determination of such an upper limit, see Ref. [521 . Fluxes of this size are required in the context 
of the Z-burst scenario for the highest energy cosmic rays. 
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Figure 3: Reach in the black hole production parameters for 5 = 6 extra dimensions, for contained 
events in an under-ice detector at a depth of 2 km and with an 1 km^ ficudial volume (1 yr = 
3.16 • 10^ s). The shaded solid lines labeled "FE" indicates the constraint |16] arising from the 
non-observation of horizontal showers by the Fly's Eye collaboration jTH]- The shaded dotted, 
My) = iWbh"^, and shaded solid, = (1/5) M^™, lines give a rough indication of the boundary 
of applicability of the semiclassical picture 0. Top: Exploiting the cosmogenic neutrino flux from 
Ref. |42i (cf. Fig. Ej). Bottom: Exploiting the upper limit on the neutrino flux from "hidden" 
hadronic astrophysical sources from Ref. |44j (cf. Fig. |2)). 
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Table 1: Event rates (1 yr = 3.16 ■ 10^ s) of contained reactions and through-going muons, 
respectively, for Standard Model like cross sections for the case of the cosmogenic neutrino flux 
from Ref. assuming a maximum energy of -Emax = 3 • 10^^ eV for the ultrahigh energy cosmic 
rays, and the case of the upper limit on the neutrino flux from "hidden" hadronic astrophysical 
sources from Ref. [Hj, respectively. An under-ice detector at a depth of 2 km has been assumed. 
For the contained events, we have considered, as in Eq. (|7j), charged current reactions from 
all flavors, but have neglected neutral current reactions, because the fraction of visible energy, 
{y) ~ 0.2, as well as the cross section, is lower for the latter [SI]. For the through-going events, 
we have considered charged current Uf^N and u^N interactions and an equipartion of neutrino 
flavors in the neutrino flux. 



LHC starts operating and Auger has taken data for already a few years. For the under- water /ice 
neutrino telescopes considered so far, the best perspective for detection of black hole events on 
the basis of the contained events might thus be assigned to RICE |2ZI, a currently operating 
radio- Cherenkov telescope with an effective volume of about 1 km^ for 10^ GeV electromagnetic 
cascades 



4- The ability to detect muons from distant neutrino reactions increases an under-ice/ water 
detector's effective neutrino target volume dramatically and is the premise upon which such 
detectors can act as neutrino telescopes. Let us therefore, in addition to the rate © of fully 
contained black hole events, estimate the rate of through-going muons at a neutrino telescope. 
We define through-going muon events to be events with a vertex located outside of the detector 
and which then pass through the detector. The expected rate for such events is then proportional 
to the area of the detector. The muon rate above some threshold i?th per unit area A and solid 
angle fl can be estimated by 



d^iV^h 
dAdtdQ 



°° 1 r'^ r 
J2 / dE,F,{E,) dn, 

TTLn JEa, J 



n. 



dn,, k ! 



X 



(9) 



X 



x{e) 



dX'exp -al%{E,){Xi9)-X')/mM [1 - [l ~ pf{E,, Et^, X' 



In words, the rate of through-going muon events from a black hole reaction is obtained by summing 
over different muon multiplicities k, weighted with the Poisson probability of mean that k 
muons emerge from black hole decay and that at least one can be observed in the detector. The 
quantity p^, which may be interpreted as the detection probability of a "typical" muon produced 
in a black hole event, is defined as (cf. also Ref. pTHj ) 



pfiE,,E,^,X')= r dE^ 



1 dn^^ 

^ p,{E„E,^,X') 



' < dE, 



(10) 



with p^{E^, Eth_, X') denoting the probability that a muon of definite energy E^ traverses an 
amount X' of material with Efj_ > Eth- 
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Figure 4: Reach in the black hole production parameters for (5 = 6 extra dimensions, for through- 
going muons in an under-ice detector at a depth of 2 km and with an 1 km^ effective area 
(1 yr = 3.16- 10^ s). The shaded solid lines labeled "FE" indicates the constraint arising from 
the non-observation of horizontal showers by the Fly's Eye collaboration jTH|. The shaded dotted, 
Mu = M^"^, and shaded solid, = (1/5) M^™, lines give a rough indication of the boundary 
of applicability of the semiclassical picture 0. Top: Exploiting the cosmogenic neutrino flux from 
Ref. |42i (cf. Fig. Ej). Bottom: Exploiting the upper limit on the neutrino flux from "hidden" 
hadronic astrophysical sources from Ref. |44j (cf. Fig. |2)). 
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If one neglects the stochastic effects of range stragghng which may become important for muon 
energies close to the detection threshold, i.e. within the approximation of a deterministic energy- 
range relationship, one has the step function 



p^{E^,E,i„X')=eiE^-E'(E,^,X')) (11) 



where E^,, is the solution of 



x'=r''^^^— (12) 

Je,^ a{E) + b{E) E' ^ ' 

The ionization and radiation loss parameters, a{E) and b{E), vary only slowly with energy and 
take the asymptotic values a ~ 3 ■ 10^^ GeV cm^/g and 6^4- 10~^ cm^/g jSHj- Choosing these 
constant values, we get for the average muon range. 



^th+ 1-e 



-bX' 



e"-^ . (13) 



Since for our considerations -Eth ~ 10^ GeV, we can safely neglect the contribution proportional 
to a in the following. 

For the laboratory energy distribution of the muons produced in a black hole event, entering the 
probability ^U^. we have taken for simplicit}0 (cf. also Ref. jH^ ) 



n 



l^dE^ 2E, \{n) % 



with (n) being the average number (0) of particles produced in black hole decay (cf. Fig. The 
distribution (fT^ is obtained by boosting a delta-function energy spectrum. 



dEf,"^ 



oc5(i?r-Mbh/(T2)), (15) 



in the cm frame of the black hole, with isotropic direction of the decay particles, with a Lorentz 
factor 7 = E^/M^Yl to the laboratory frame. 

Finally, we have exploited the ansatz (similar to Eq. p5|) ). 

= a,^]^ 6 (n, - {n,)) , with {n,) ^ {n)/30 , (16) 

to perform also the integration over the multiplicity in our rate estimate Q. The above 
relation for (n^) is based on the assumption that a black hole does not discriminate between 
Standard Model particles in its evaporation spectrum j;8j. 

The projected reach in the black hole production parameter space for 5 = 6 extra dimensions 
for through-going muons in an under-ice detector at a depth of 2 km and with an effective area 

^Alternatively, one may use a boosted thermal distribution, ^ Se" /2-f ^E{E /TyiY / {e^ /'^^ + 1), which, however, 
needs an additional numerical integration. 
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of 1 km^ is shown in Fig. |3J for two of the discussed ultrahigh energy cosmic neutrino fluxes. 
The corresponding expected Standard Model background rates for through-going muons due to 
charged current interactions is summarized in Tabled (see also Ref. 57 for related calculations). 
An additional source of background which is not related to an astrophysical neutrino flux is the 
muon flux arising from charm production by cosmic rays interacting in the atmosphere. Using 
two available models |58[ IKH] we have estimated a background of 0.33 ^1.5 through-going muon 
events from atmospheric charm production per year. To further reduce this background, one may 
introduce a cut in the zenith angle, e.g. 6* > 60°. 

5. We have considered the sensitivity of neutrino telescopes to black hole production in the 
context of extra dimension scenarios with TeV-scale gravity and compared it to the one of air 
shower facilities and the LHC For concreteness, we have assumed a detector with 1 km^ area and 
1 km^ volume, respectively, located at a vertical depth of 2 km in ice. 

We have seen that the effective volumina for contained events, V ~ 0.001 ^ 0.01 km^, available 
from the currently operating AMANDA and Baikal detectors are too small to yield a large enough 
rate to be competitive with the already existing constraints on black hole production from Fly's 
Eye and AGASA and with the projected ones from Auger ^21 ^1 El- The planned IceCube 
detector with a fiducial volume of about 1 km^ would have an interesting perspective. However, 
with the start of construction of IceCube in 2004 and full completion only after the physics start 
at the LHC, the discovery potential for black hole production is strongly reduced. Instead, the 
radio- Cherenkov telescope RICE [2Z|, with an effective volume of about 0.1 (1) km^ for 10 GeV 
electromagnetic cascades [53], could set sensible constraints using already available data. 

The prospects for a search for black hole signatures using through-going muons are also good. For 
a detection area of 1 km^ we have shown that in case the neutrino flux is just on the level of the 
cosmogenic one, only a few (<1) contained or through-going background events from Standard 
Model charged current interactions are expected per year. The background of prompt muons 
from atmospheric charm was estimated to be of similar magnitude and can be further reduced 
by application of a cut on the zenith angle. It was shown in Ref. [^0]; that AMANDA achieves 
an effective area of about 0.3 km^ for down-going muons above an energy of 10^ GeV and that 
suppression of the atmospheric muon background is possible. With 5 years of AMANDA data 
available and an assumed duty-factor of 2/3 this leads to an exposure of about 1 km^-yr. There- 
fore, strong constraints on the black hole production parameters are expected from AMANDA 
if the currently available data show no through-going muon^ above lO'' GeV. In this case, the 
sensitivity of AMANDA - approximately the 3 events/yr contour in Fig. 0] (top) for a 90 % confi- 
dence level limit - compares favorably with the one of Fly's Eye, AGASA, and Auger [T2 t ll6 | IT7j. 
And even at the kinematic limit of the LHC, MJJJJ'^ 10 TeV (cf. Ref. p£,), one is still left with 
a sensitivity on Md of about 1.3 TeV. 

In the case that the ultrahigh energy cosmic neutrino flux considerably exceeds the cosmogenic 
one, discrimination between Standard Model background and black hole events becomes crucial. 
This might not be achieved by observation of through-going events alone, but with the addi- 
tional information of contained events. For example, the spectral information from contained 

^It is noteworthy that in this case AMANDA can put also strong constraints on sphaleron production param- 
eters PS]; since the phenomenology of black hole and sphaleron production are quite similar. 
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and through-going events should be remarkably different, as the through-going muons from black 
hole decays are produced with considerably lower energy. Another discriminator may be based 
on the different hadronic components of Standard Model processes and black hole events. In a 
typical Standard Model charged current interaction one has an average hadronic energy fraction 
of (y) ~ 0.2, whereas black hole events have typically a large hadronic component (5 : 1) [8^ corre- 
sponding to a {y) ~ 0.8. In particular radio- Cherenkov detectors such as RICE have the potential 
to discriminate between electromagnetic and hadronic showers by measuring the characteristically 
different angular radiation patterns [HTj . 

To summarize, we have shown that large scale under- water/ice neutrino detectors have a sub- 
stantial potential to discover or rule out black hole production which is comparable to that of 
air shower facilities. In the optimistic case that a neutrino flux significantly higher than the 
cosmogenic one is realized in nature, one even has discovery potential for black holes at neutrino 
telescopes beyond the reach of LHC 
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